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I review the status of non-perturbative investigations of the finite temperature electroweak phase transition by 
means of lattice simulations. 



1. Introduction 

Since it has been suggested [1] that the baryon 
asymmetry of the universe could have been gener- 
ated at the symmetry restoring [2] finite tempera- 
ture electroweak phase transition, much effort has 
been invested to study this idea in detail. If the 
above scenario is realized, the question of how the 
observed baryon asymmetry has been generated 
could be answered within the framework of the 
minimal standard model of the electroweak inter- 
actions. It is therefore a natural first approach to 
use perturbation theory for studying this ques- 
tion. Indeed, as of today results up to the 2-loop 
level are available [3-5,24]. However, it soon has 
become clear that even with the incorporation of 
special techniques to cure the infrared problems 
occuring in the symmetry restored phase, the per- 
turbative series show bad convergence properties 
at least for Higgs masses larger than about 50GeV 
[5]. As a consequence one may question the re- 
liability of the results as obtained in perturba- 
tion theory. In addition, it was found that for 
Higgs masses Mh ~ 80GeV the perturbative se- 
ries break down. Therefore only a small window 
of possible Higgs masses is left since the present 
experimental lower value is Mh > 64GeV [6]. 
This window would even be closed if one takes 
theoretical lower bounds for the Higgs mass as a 
result from vacuum instability [7]. Nevertheless, 
it would be very important to know whether per- 

* Plenary talk given at the International Symposium 
on Lattice Field Theory, 11-15 July 1995, Melbourne, 
Australia 

^e-mail: kjansen@desy.de 



turbation theory can provide a quantitative de- 
scription of the electroweak phase transition up to 
Higgs masses of about SOGeV. A positive answer 
would give strong confidence that we understand 
the properties of the electroweak phase transition. 

It is obvious that only non-perturbative meth- 
ods like numerical simulations can shed light on 
the reliability of the results as obtained in per- 
turbation theory. Since the heart of the infrared 
problems arising in perturbation theory is the 
SU(2)-Higgs part of the minimal standard model, 
one can in a first approximation neglect the U(l) 
sector and the fermions. In this setup the problem 
becomes well suited for a numerical study since 
simulations of the SU(2)-Higgs model on the lat- 
tice are straightforward. In particular one does 
not have to deal with numerical simulations of 
chiral fermions on the lattice, a problem that as 
of today has not been solved [8]. 

Out of the three conditions for baryogenesis as 
formulated by Sakharov [9], -baryon number and 
CP violation and out of equilibrium processes- 
lattice simulations can test the last condition. 
The reason is that a possible realization of non- 
equilibrium behaviour is a first order phase trans- 
ition. Lattice simulations are able to answer the 
question of the nature and -in case of a not too 
weak first order phase transition- the strength of 
the finite temperature phase transition. Whether 
the right amount of CP violation is generated will 
not be addressed here. 

Any baryon asymmetry that has been gener- 
ated at the electroweak phase transition is washed 
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out with a rate 



adopted. The standard lattice form of the SU(2)- 
Higgs model is given by the action 



(1) 



where g is the electroweak SU(2) gauge coupling, 
TP'^ is the physical temperature and v is the scalar 
vacuum expectation value to be taken at tem- 
peratures TP'^ ^ Tc with Tc the physical transi- 
tion temperature of the electroweak phase transi- 
tion. Although the coefficient in the exponent is 
not known exactly, semi-classical estimates for its 
value exist [10]. Taking these estimates and com- 
paring the rate in eq.(l) to the expansion rate 
of the universe leads to the rather conservative 
bound 



v{TP^)/TP^>l 



(2) 



in order that baryon asymmetry is not washed 
out [11]. This implies that the phase transition 
not only has to be of first order but that it has to 
be strong enough, with "strong enough" defined 
by relation (2). 

The values of the Higgs mass that have been in- 
vestigated so far in numerical simulations are less 
than Mh = 70GeV. The reason for the restric- 
tion to these small Higgs mass values is that the 
strength of the phase transition weakens rapidly 
with increasing Higgs mass. For larger Higgs 
masses it becomes very difficult to resolve the or- 
der of the phase transition. Therefore all numer- 
ical work that is presented in this review mainly 
serves the purpose to confront results obtained in 
perturbative computations with the ones deter- 
mined non-perturbatively by means of lattice sim- 
ulations. Earlier review articles concerning the 
lattice approach to the finite temperature elec- 
troweak phase transition were given by Shaposh- 
nikov [12] and by Kajantie [13]. 

2. The SU(2)-Higgs model: action and 
phase diagram 

As mentioned in the introduction, the focus of 
this review will be on the SU(2)-Higgs part of 
the minimal standard model and the fermions as 
well as the U(l) part are neglected. To perform 
computer simulations a lattice regularization is 
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Here Uxp, denotes the SU(2) gauge link variable, 
Upi is the product of four [/'s around a plaque- 
tte and Lp^ is a complex 2 Cg) 2 matrix in isospin 
space describing the Higgs scalar field and sat- 
isfying = T2(p'^T2. Often the length of the 
Higgs scalar field, p'^ = ^TiLp\.Lpx, will be used. 
The plaquette term assumes the form of the pure 
Yang Mills action in the continuum limit when 
the lattice spacing a is sent to zero. The param- 
eter [3 is related in the usual way to the SU(2) 
gauge coupling g hy j3 = A/g'^ . The terms on the 
second line are the scalar bare potential with bare 
quartic coupling A. The last term is the gauge in- 
variant kinetic part of the scalar fields. It assumes 
the continuum expression of the gauge covariant 
derivative when a ^ 0. The scalar hopping pa- 
rameter K is related to the bare mass mo by hIq = 
(1 — 8k — 2X)/k and the quartic coupling A to its 
continuum counterpart by Aq = A/4k^. Tuning k 
is essentially equivalent to tuning the bare mass 
parameter. In the following, the lattice spacing 
a will be set to one unless otherwise stated. To 
realize a finite temperature the time extension Li 
of the lattice with volume V = LiL^LyL;, is peri- 
odized and kept much smaller than the the spatial 
directions L^^y^z which are chosen large enough to 
mimic infinite volume lattices. 

At zero temperature considerable effort has 
been expended on the SU(2)-Higgs model in 
numerical lattice simulations [14]. The situa- 
tion at a non-zero temperature is quite differ- 
ent and much less is known there. The pioneer- 
ing works of [15] and [16] concentrated on large 
Higgs masses, that is on large values of the bare 
quartic coupling. The symmetry restoration pic- 
ture could be verified in these works. It was also 
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shown that the lattice data are consistent with an 
analytical crossover behaviour of the electroweak 
phase transition at large values of A. 

The results of these studies give the following 
schematic picture of the phase diagram in the 3- 
dimensional parameter space given by the cou- 
plings /3, K and the temperature aT'^^ = T = 
l/Li depicted in fig. 1 for a fixed A cx 0(1). 
The phase diagram should be understood as be- 
ing obtained on lattices with infinite extensions in 
the spatial directions and finite time extensions 
such that the highest temperature is T = 1 cor- 
responding to _L( = 1. 




Figure 1. Schematic finite temperature phase di- 
agram of the SU(2) Higgs model. 

For K = we are left with the pure SU(2) gauge 
theory. This model exhibits at any T > the fi- 
nite temperature deconfinement phase transition 
at a critical value /3c which shifts to larger /3 val- 
ues when T is decreased. Switching on the scalar 
fields by taking a nonzero value of k this phase 
transition is stable and extents into the cube. 
For K large enough the transition may change to 
an analytical crossover behaviour as suggested in 
[16]. For l3 = CO the gauge fields are pure gauge 
and may be set to unity. The resulting model 



is the theory with global SU(2)i(g)SU(2)ij or 
0(4) symmetry. Keeping T > fixed while vary- 
ing K the electroweak finite temperature phase 
transition is passed at a critical Kg. This phase 
transition separates the 0(4) symmetric from the 
symmetry broken phase with an residual 0(3) 
symmetry. Investigations of this model revealed 
only second order phase transitions [17-19]. The 
same is true if fermions with vector-like Yukawa 
interactions are added [20]. Releasing the tem- 
perature the finite temperature phase transition 
becomes a phase transition line. Choosing also 
l3 < CO this develops into a transition surface. 
For small /3 one may, similar to the zero tem- 
perature case, expect a region in the parame- 
ter space where observables behave analytically 
and no phase transition takes place. For A and 
temperature large enough, also the electroweak 
phase transition may change into a crossover phe- 
nomenon. In summary there are three regions in 
the phase diagram: 

• confinement region for /3 < /3c, k ^ 1 

• Higgs region for /3 ^ 1 and k > Kg 

• deconfinement region for [3 > [3^, k < 

To make contact to continuum results like the 
ones obtained in perturbation theory, the lattice 
study has to control finite size effects and finite a 
effects. Since a physical temperature is related to 
the time extent of the lattice by T'^^ = 1/Lta, the 
continuum limit a ^ corresponds to send the 
time extent Lt ^ co. If one considers the case of 
finite temperature QCD [29] the necessary finite 
time extent of the lattice as used in the numerical 
simulations produce a severe limitation showing 
up in strong finite a effects. This implies that 
large values of Lt have to be taken. Although it 
will be shown below that the situation in the elec- 
troweak theory is different, this danger motivates 
the use of the dimensional reduction program that 
will be described in the next section. 

Since the different parts of the phase diagram 
may be analytically connected it would be more 
appropriate to call them regions. However, in the 
literature -and sometimes in this article- they are 
also denoted as Higgs, confinement and deconfine- 
ment (or symmetric) phases. It is important to 
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note that the finite temperature phase transition 
occurs at /3 = 8 — 10 if one takes the weak SU(2) 
gauge coupling as input. Therefore the transition 
occurs between the Higgs region and the decon- 
finement region. When T is decreased to reach 
the continuum limit, while keeping the physical 
situation fixed, the finite temperature transition 
shifts to larger /3- values. Thus the electroweak 
finite temperature phase transition always takes 
place between the Higgs and the deconfinement 
regions and does not interfere with the remnant 
of the pure gauge transition. 



3. Dimensional reduction 

In the euclidean path integral formulation of 
quantum field theories a finite temperature is in- 
troduced by periodizing in time. The Matsubara 
frequencies generated by this procedure can be 
separated into two kinds: The static Matsubara 
frequencies corresponding to massless zero modes 
and the non-static ones corresponding to massive 
modes. Since it is the zero Matsubara frequen- 
cies that lead to infrared problems in perturba- 
tion theory, it is a natural idea to integrate out 
all non-static modes. In this way the original 4- 
dimensional theory at finite temperature will be 
reduced to an effective 3-dimensional model. Of 
course, this will not cure the infrared problems 
associated with the static Matsubara frequencies. 
The idea in the approach of dimensional reduc- 
tion is now to perform a simulation of the ef- 
fective 3-dimensional model and treat the non- 
perturbative content numerically. The program 
as sketched above has been thoroughly discussed 
and carried out in [21-25] (for latest results, see 
[27,26]). 

In practice the integration can not be per- 
formed exactly but the reduction is done in per- 
turbation theory and one may wonder when this 
procedure is valid. In [25] this question has been 
investigated in detail. The reduction itself is jus- 
tified when perturbation theory ai zero tempera- 
ture is valid. More specifically, the 4-dimensional 
finite temperature effective potential, V^^^^ , is 
approximated by its 3-dimensional counterpart. 



V^J^j: , through 



o 



(4) 



In (4) m stands for a relevant mass scale in the 
zero temperature 4-dimensional theory. Since 
there all masses are proportional to the couplings 
A and g, dimensional reduction is valid when 
the 4-dimensional (renormalized) zero tempera- 
ture couplings are small, i.e. when zero tempera- 
ture perturbation theory works. 

However, this does not imply high temperature 
perturbation theory to work. For this an addi- 
tional condition has to be satisfied, namely that 
the ratio gs/niT = g'^T'^^ / rriT ^ 1 with g^ the 3- 
dimensional gauge coupling and my a finite tem- 
perature mass scale. Performing the reduction 
perturbatively to 1-loop order one finds the fol- 
lowing effective lattice action: 
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+ ^Tr^ _ ip'^{x)Ui{x)ip{x + i)) 

1 - 2i3r - 3i3h u \ r \ 
+ 2^ 2 ^T:'pHx)ip{x) 

X 

-^E^TrAg(a.)iTr^t(^)^(^) 

X 

+ /5«EQTr^^W^w) • (5) 

In (5) the time component Aq of the gauge field 
(in continuum notation) has been left in. In prin- 
ciple it can also be integrated out. The constant 
S is given by S = 0.252731. Summations are 
to be taken over 3-dimensional lattices. The ac- 
tion contains the 4-dimensional gauge coupling g 



5 



and three additional bare couplings [Sq , Ph and 
Pr. They have a well defined relation to the 4- 
dimensional bare couplings. The coupling [Sq is 
given by 

4 1 

= —77;^ • (6) 

Eq.(6) implies that for fixed and physical tem- 
perature TP'^ tuning /3g ^ 00 corresponds to the 
continuum limit a ^ 0. I3r is given by 



8M^ f^G 



(7) 



with Mh and Mw the zero temperature 4- 
dimensional Higgs and W-mass, respectively. The 
parameter i3h is left free and has to be tuned to 
the transition point. 

As input 4-dimensional parameters, the gauge 
coupling g = 2/3, the W-boson mass Mw = 
80.6GeV and as a free parameter the Higgs mass 
Mh, is chosen. Note that the Higgs mass is a 
renormalized MS mass to be taken at a scale 
yjipft physical value might be different. 

The numerical computations obey the follow- 
ing strategy. One first fixes i3g which fixes the 
lattice spacing a. Then for a given lattice size i3h 
and with the help of relation (7) I3r are tuned 
to the transition point. The finite lattice tran- 
sition points are determined by various methods 
(see [27] for a description of the methods applied). 
Although these methods may yield different /Sfj 
on the finite lattice, the infinite volume limit ex- 
trapolation should give the same /Sfj. Therefore 
a careful finite size scaling is performed to obtain 
l3'fj at infinite 3- volume Vsd- This procedure is 
repeated for larger and larger i3g and hence for 
smaller and smaller a values. 

For every fixed set of parameters and with the 
non-perturbatively determined value of jS'fj phys- 
ical quantities can be computed. One example is 
the ratio Mh/Tc which is given by 
2 r 
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Ml 
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(8) 
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In the above formula which is calculated in per- 
turbation theory every quantity beside the con- 
stants rj, fj and i) is given. The constants can be 
fixed by relating the lattice scheme to a contin- 
uum scheme chosen to be MS. In earlier works 
the constants have been evaluated by computing 
condensates like < ^Tr(pl,(px > in the MS scheme 
and measuring it directly on the lattice. In the 
relation of the condensates computed in both 
schemes the unknown constants appear. Measur- 
ing < ^Tifl.Lpx > at various coupling parameter 
values allows therefore for a fit from which the 
constants can be extracted. 

Recently, Laine [28] performed a full 2-loop per- 
turbative computation relating the lattice to the 
MS scheme. This allows for a complete 2-loop 
perturbative computation of the numerical values 
for the above mentioned rj constants. Comparing 
these 2-loop results with the earlier Monte Carlo 
results for the constants a satisfactory agreement 
was found. Taking eq.(8), an infinite volume -but 
still non zero lattice spacing- value for Mh /Tc 
can be given. By repeating the numerical sim- 
ulation at larger values of i3g , and therefore at 
smaller lattice spacings, finally the continuum 
limit extrapolation can be achieved. 

The procedure is illustrated in fig. 2 and fig. 3. 
Fig. 2 shows the infinite volume extrapolation of 
I3'jj at fixed Mh = 35GeV and fixed /3g = 8. The 
figure nicely demonstrates the convergence to the 
same infinite volume /Sfj from all methods used. 
Fig. 3 shows the extrapolation to the continuum 
a = 0. Note the scale in the plot which indicates 
the very precise determination of Tg. 

Progress since the Bielefeld conference last year 
in the 3-dimensional reduced model consists in 
two things. First, there is now a full 2-loop ana- 
lytical expression of relating lattice data to their 
continuum counterparts in the MS scheme. Sec- 
ond, new simulation results have been obtained 
at Mh = 60 and Mh = 70 GeV. Difficulties 
to resolve the order of the phase transition at 
Mh = 80GeV were reported at this conference 
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Figure 2. Infinite volume extrapolation for jS'fj 
from four different methods as explained in [27]. 

[27]. In table 1 a comparison is made between the 
transition temperature as obtained from the sim- 
ulations in the 3-dimensional reduced model and 
perturbation theory done in the effective model in 
3 dimensions. As a second quantity for compar- 
ison the ratio v(Tc)/Tc is taken. Note that the 
data in table 1 can not be directly related to 4 
dimensional values. 

Table 1 demonstrates the high precision that 
can be obtained in simulations of the SU(2)-IIiggs 
model. Differences between the perturbatively 
obtained results and from the simulations can be 
seen. But these discrepancies are very small indi- 
cating also only small non-perturbative effects. 

3.1. Critique on dimensional reduction 

Recently, Jakovac [30] computed the 3- dimen- 
sional reduced theory up to 2-loop order in a 
scalar field theory. The effective model in three 
dimensions is super-renormalizable. All counter 
terms that are generated in the reduction step 
have to be cancelled by terms appearing in this 
super-renormalizable theory. Although this cri- 
terion works perfectly at the 1-loop level, a clear 
mismatch is encountered at 2-loop order. The 
counter terms generated in the reduction process 
can no longer be cancelled by the ones in the 3- 
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Figure 3. Continuum limit extrapolation of the 
critical temperature Tg given in GeV. 

dimensional theory, see also [31]. 

Cancellation can only be achieved again if non- 
local terms in the effective action are allowed. 
These non-local terms are expected on quite gen- 
eral grounds to appear on the scale (2'itT)~^ [32]. 
Splitting the low frequency from the high fre- 
quency modes and treating the low frequency 
modes as the block averages [33] one obtains an 
expression for the finite temperature propagator 



TT(r,t) = T 



(27r)3 



-ipr — iLOni 



1 



(9) 



which leads to exponential decay modes with de- 
cay length (27rT)~^. To circumvent this problem, 
in [30] it was suggested to integrate out all modes 
with p^ T and to obtain in this way an effective 
cut-off theory. A different approach was followed 
in [32]. There it was suggested to stay on the 
lattice and to perform block-spin transformations 
-either analytically or numerically- to separate 
the low from the high frequency parts. For finite 
temperature systems the block size in the tem- 
poral direction should be taken to be 1/T. The 
remaining (perfect) effective block spin action is 
then purely 3-dimensional. Results for the elec- 
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Table 1 



Comparison of the perturbatively [pert] and numerically [latt] obtained results in the 3-dimensional 
reduced model. Mh and Tg are given in GeV. 



Mh 


Tc [pert] 


Tc [latt] 


v{Tc)/Tc [pert] 


v{Tc)/Tc [latt] 


35 


93.3 


92.64(7) 


1.75 


1.86(3) 


60 


140.1 


138.20(4) 


0.68 


0.691(7) 


70 


157.0 


154.5(1) 


0.55 


0.57(2) 



troweak phase transition emerging from this at- 
tractive concept have still to be awaited for. In a 
similar spirit is the average action approach in the 
continuum [34]. In the approach of dimensional 
reduction the problem of the non-local terms is 
treated by matching Green's function of the 3- 
dimensional and 4-dimensional theories [24,26]. 

4. Four dimensions 

The conceptually cleanest and clearest way to 
treat the finite temperature SU(2) Higgs model is, 
of course, a direct simulation of the four dimen- 
sional theory at finite temperature. However, as 
Kanaya's talk at this conference again confirmed, 
the lesson from finite temperature QCD simula- 
tions is that one might have to expect sizable scal- 
ing violations. This would result in large values 
of Li to suppress finite a effects which in turn ne- 
cessitates large lattices since one wants to keep 
Li Lg. This danger motivated the dimensional 
reduction program. It seems, however, that we 
are fortunate in the SU(2)-IIiggs model. As will 
be shown below, scaling violations for the quan- 
tities of interest are surprisingly small and safe 
extrapolations to continuum physics are possible. 

Lattice simulations of the 4-dimensional theory 
have been performed by two groups. In [35] the 
finite temperature simulations at small A were ini- 
tiated. Newer results of this group were presented 
last year at the Bielefeld lattice conference [36,37]. 
Results for a determination of static potentials 
can be found in [38]. A high statistics simulation 
of the 4-dimensional SU(2) Higgs model was per- 
formed in [39-42]. Due to the use of the Alenia 
Quadrics (APE) massively parallel machine and 
of algorithms specifically designed for the SU(2)- 
Higgs model [43,36] a high statistics and precise 
data could be obtained. It should be noted that 



the realization described in [36] of the method 
introduced in [43] gives the best results for the 
autocorrelation time. Fortunately enough, both 
groups performed simulations at the same param- 
eter values. The results they obtain for various 
quantities are in good agreement. In the follow- 
ing I will therefore concentrate on the methods 
and results of only one group. 

The basic goal of all the simulations is to make 
contact to continuum physics. To reach the con- 
tinuum limit one has to run along renormalized 
trajectories also called lines of constant physics. 
Of course, the a = limit can never be reached 
in the SU(2)-Higgs model because of triviality. 
What we mean by continuum limit therefore is 
to approach the scaling region of the SU(2) Higgs 
theory where cut-off effects are negligible. There 
the approximation of the continuum cut-off the- 
ory is expected to be very good in the low energy 
regime. The lines of constant physics are given by 
dimensionless ratios of physical quantities or the 
dimensionless coupling constants gn, the renor- 
malized gauge coupling and \r = M^/Sv^, the 
renormalized quartic coupling. It is natural to fol- 
low perturbation theory and perform the renor- 
malization at zero temperature. 

In the 4-dimensional approach to the finite 
temperature electroweak phase transition in a 
first step the bare parameters of the action (3) 
are fixed such that the desired physical situation 
is realized. As physical input ones takes the W- 
boson mass Mw = 80GeV and the gauge cou- 
pling /3 = 8 which amounts to = 0.5. Fixing 
the Higgs mass basically amounts to fixing the 
bare quartic coupling A. Choosing A = 0.0001 the 
Higgs mass will be around 20GeV and A = 0.0005 
corresponds to Mh ~ 50GeV. The physical tem- 
perature TP'^ = T/a is given by the temporal ex- 
tent Li of the lattice. The phase transition itself 
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is then found by tuning the scalar hopping param- 
eter K to its transition value Kg (see fig. 1). At Kg 
one determines the order of the phase transition. 
If it is first order, quantities like surface tension, 
latent heat, jump of the order parameter etc. are 
computed. To get the values for the renormal- 
ized zero temperature quantities additional sim- 
ulations at K = Kc are needed but now at zero 
temperature which amounts to choose a hyper- 
cubic symmetric lattice which is large enough to 
suppress finite size effects. Note that k = Kg at 
zero temperature belongs always to the symmetry 
broken Higgs phase, see fig. 1. 

The approach to the continuum limit a ^ is 
realized by increasing Li while staying on the lines 
of constant physics. To compensate the change of 
a the bare parameters /3 and A have to be changed 
accordingly. To estimate the values of these pa- 
rameters when going from Li to _L( -|- 1 the 1-loop 
continuum renormalization group equations are 
used 



dr 



dXo(T) 
dr 



43 
"48^ 



167r2 



0{\l,\ig\ 
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(10) 



with r = hi{aMw)~^ ■ 

From the above discussion it should become 
clear that a first important ingredient in the study 
of the electroweak phase transition is the determi- 
nation of the transition k^. For this three meth- 
ods have been used. It is noteworthy that in the 
4-dimensional simulations for the determination 
of all relevant quantities at least two independent 
methods are used, providing very useful cross 
checks. The most precise method to determine 
Kc is the multicanonical simulation technique [44] 
which gives precision data for the distribution of 
suitable order parameters. Applying the equal 
height or equal weight criterion to these distri- 
butions allows to give up to seven digits [42]. 
The presumably most practical method is the so 
called 2 — K method [40,42]. Here one chooses 
elongated lattices and sets one half of the lattice 
to Ki belonging to the symmetric and the other 



half to K2 belonging to the Higgs region of the 
model. In this way an interface is generated be- 
tween the two halfs. Bringing both k values closer 
to Kc, eventually the interface will break down. 
The corresponding ki and K2 give a lower and 
upper bound for the transition k^. This method 
works to a good and satisfactory precision and is 
mostly used for a first estimate of k^. 

As an analytical tool the effective potential can 
be used. In [40] the gauge invariant effective po- 
tential [45,46] has been used. It was recently 
found that at least in 1-loop order the effective 
potential in Landau gauge is equivalent to the 
gauge invariant one if one chooses a particular 
mass resummation in the Higgs region by setting 
the Goldstone boson mass to zero [5]. One ob- 
tains for the broken phase 



^l — loop — ^tree 

d'^k fg 



(11) 



where the masses are related to the parameters 
in the lattice action (3) and the translational in- 
variant saddle point solution [5] by 



,2 „2 



(12) 



and the momenta in the lattice integrals are k'^ = 
^^[2 — 2cos(A;^)]. The tree level potential reads 

Vireeip'') = (1 " ^n)p'' + A(p2 - 1)2. (I3) 



The Landau gauge effective potential allows 
for a standard renormalization. One defines the 
renormalized mass M'^/Zr = d'^V/dp'^\^_^ 
and the renormalized coupling by Xr/Z"^ = 
d'^V/dp'^\^_^ . Here Zr is the wavefunction 
renormalization constant [47]. The renormalized 
vacuum expectation value is vr = vJ^JZr. For 
the value of the gauge coupling used in evalu- 
ating the effective potential the mean field im- 
provement is taken, g'^ g'^ j < Up > [48] with 
< Up > the averaged measured plaquette value. 
When below lattice simulation results are com- 
pared to lattice perturbation theory, the effective 
potential (11) and the above definitions of renor- 
malized quantities are used. Note that in the ef- 
fective potential approach for the renormalization 
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step at zero temperature the value of Kg is taken 
as determined from the effective potential itself. 

All three methods mentioned above give a first 
order phase transition. The most clean signal 
comes from the distributions of suitable observ- 
ables using multicanonical simulation techniques. 
They show very pronounced two peak signals with 
about 10 orders of magnitude difference between 
the maximum and the minimum of the distribu- 
tion [40]. After having identified the first order 
nature, as a second step computations of quanti- 
ties characterizing the strength of the phase tran- 
sition are necessary. A first such quantity is the 
latent heat. 

4.1. Latent heat 

Its calculation may proceed in the following 
way [5]. At the phase transition the free energy 
W is continuous and 



W,(T,.J(t)) = Wi,(T,J(T)) 



(14) 



determines the phase boundary. Here Ws(b) is 
the free energy of the symmetric (broken) state. 
The external source J is thought of being cou- 
pled gauge invariantly to the composite operator 
Tr(p^(p. Taking the total variation in eq.(14) one 
obtains 



d.J(T) 
dT 



(15) 



which gives the latent heat AQ through -^{Wg 



m) 



-AQ/T and the order parameter jump 



Ap"^ via --§j(Ws - Wi) = Ap'\ As a result the 
well known Clausius-Clapeyron equation 

d.] 



AQ = Ap'T 



dT 



(16) 



is obtained. 

For the model under consideration for dimen- 



sional reasons we have -\- J{T) 
And therefore 
dj 



T 



dT 



2mi 



C{g\\)T-' 



(17) 



j=o 
and finally 

AQ = -M^Ap^l 

with Mh 



-0{g\X)) (18) 

2mg. This relation connects the 
order parameter jump to the latent heat. A sim- 
ilar expression was obtained and used also in the 



3-dimensional reduced model. There is a second 
more complicated expression for the latent heat 
that can be obtained by differentiating the action 
density with respect to the lattice spacing. This 
second definition is explained and used in [40]. 



<1 



0.1 




20 40 
Mjj [GeV] 



60 



Figure 4. Latent heat as function of the zero 
temperature Higgs mass Mh- The full symbols 
are the definition using the derivative of the ac- 
tion density with respect to the lattice spacing 
[40]. The open symbols use the definition by the 
Clausius-Clapeyron equation (16). The diamond 
shaped symbols are from lattice perturbation the- 
ory. Squared symbols denote Lt = triangles 
Lt = 2. 

In Fig. 4 a comparison is made for both def- 
initions of the latent heat. In the same figure 
results from lattice perturbation theory are plot- 
ted. There are several important features worth 
mentioning in the graph. First one notices that 
for Mh ~ 20GeV both definitions and the re- 
sults from lattice perturbation theory agree very 
well. At the higher Higgs mass again the two def- 
initions coincide but lattice perturbation theory 
does not describe the data as well any more. Sec- 
ond, there is a rapid decrease of the values for the 
latent heat when the Higgs mass is increased indi- 
cating that the phase transition becomes rapidly 
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weaker. Finally, there is almost no difference be- 
tween the Li = 2 and the Li = 3 data. The scal- 
ing violations seem to be very small which will be 
confirmed by other quantities later again. 

4.2. Surface tension 

Another important quantity characterizing a 
first order phase transition is the surface tension. 
Three independent methods to compute the sur- 
face tension will be shortly described and com- 
pared in this section. The first utilizes distribu- 
tions of suitable observables like the action den- 
sity. The distributions will develop a two peak 
structure at a first order phase transition. Of 
great help to obtain the distributions is the mul- 
ticanonical simulation technique with reweighing 
[44,49]. If by Pmax and Pmin the maximum and 
the minimum of the distributions are denoted, 
then a finite volume estimate of the surface ten- 



1 



(TV 



P 

1 fy 



P 

1 fy 



(19) 



2LxLyLi 

Taking lattices elongated in one direction, chosen 
to be and L = L^ = Ly , the infinite volume 
value for the surface tension is given by [50] 



(TV 



^{c+-\ogL, 



\^ogL) (20) 



where c is an unknown constant that can be fitted 
by choosing different size lattices. 

As a second method the above mentioned 2 — k 
method [51] can be used. For the SU(2)-Higgs 
model it amounts to measure the link expectation 
values L^ =< ^Tiipl.Ux^fj,'-px+)i > in the two halfs 
of the lattice with different ki < Kg and K2 > 
Kg. Choosing again in the z-direction elongated 
lattices one obtains [40] 

a = L,\K,,2->^c\{Ll-Ll) . (21) 

A third method is measuring appropriate corre- 
lation functions [52-54,36,42]. When the systems 
tunnels between the two states at the first or- 
der phase transition, the lowest excitation of the 
transfer matrix spectrum is the so called tunnel 
energy Eg- It is directly connected to the surface 
tension by its exponential finite size effect 



This method was recently applied to the SU(2)- 
Higgs model in [36,42]. 



Table 2 

Comparison of the surface tension a as obtained 
from the three methods described in the text. All 
results are obtained for a _L( = 2 lattice. 



Mh [GeV] 








2 - K 


distribution 


tunneling 


18 


0.84(16) 


0.83(4) 




35 


0.065(10) 




0.053(5) 


49 


0.008(2) 







Eo = C exp {-LiL^Lya} 



(22) 



Table 2 shows that the different methods give 
compatible results. Since these methods are com- 
pletely independent from each other, the agree- 
ment strengthens the trust in the numbers ob- 
tained from the numerical simulations. Note that 
also for the surface tension a rapid decrease of 
the strength of the phase transition is observed. 
Taking numbers for lattices with larger Li al- 
most no scaling violation effects can be seen [40]. 
It should be noted that for Mh £ 35GeV, cr/Tj* 
agrees amazingly well with perturbatively com- 
puted values [39]. 

5. Comparison to Perturbation theory 

In the previous two sections details of the nu- 
merical simulations in the 3-dimensional reduced 
and in the 4-dimensional theory were given. For 
both approaches it was outlined how the contin- 
uum limit is reached. As fig. 3 indicates, in the 
3-dimensional reduced theory a linear extrapo- 
lation of physical observables in i3g and hence 
in the lattice spacing a lead to continuum val- 
ues. Quadratic forms of the extrapolation do not 
change the results much. In the 4-dimensional 
model the smallness of the scaling violations, 
which are compatible to the size of the statistical 
errors, lead to trustworthy continuum estimates. 
Since also finite size effects are under control in 
both approaches, the results are comparable to 
continuum computations like 2-loop perturbation 
theory. 
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Before a real comparison can be made, the cou- 
pling parameters used in the different approaches 
should be matched. As a reference point, the 
parameter values in the 4-dimensional theory, 
Mw = 80GeV and gj^ ^ 0.57 [40], are taken. The 
results from perturbation theory can then be re- 
lated to the simulation results by choosing these 
values as input for the 2-loop formulae. It is im- 
portant to choose also a renormalization scheme 
which is close to the lattice scheme for the com- 
parison [3]. The renormalization scheme depen- 
dence can have substantial effects. To relate the 
data from the 3-dimensional simulations to the 4- 
dimensional ones, the best the author found was 
an extrapolation for physical quantities to the pa- 
rameters of the 4-dimensional simulation using 2- 
loop formulae. The physical values for the SU(2)- 
Higgs model in the 3-dimensional reduced model 
are given in [26]. 

A first important quantity is, of course, the 
transition temperature itself. The simulation 
data for the quantity Tc/Mh are very well de- 
scribed by lattice perturbation theory [40]. It is 
therefore justified to perform the continuum ex- 
trapolation _L( ^ cxD by using lattice perturbation 
theory. To get some error estimate, the error from 
the Monte Carlo data is taken for the extrapo- 
lated value. The situation is shown in fig. 5 ^. The 
two solid curves are the continuum 2-loop pertur- 
bation theory results [5], corresponding to the two 
renormalized gauge couplings found numerically 
at the corresponding two Higgs mass values. In 
fig. 5 (and fig. 6) the error of the results from 
perturbation theory coming from higher orders is 
left out. An estimate of this error can be given 
by taking the difference between the 1-loop and 2- 
loop results (see figs. 4 and 5 in [5]). It turns out 
that the discrepancy between the numerical sim- 
ulation and the perturbatively computed results 
is of the same order as this error for the physi- 
cal quantities considered in this section. Figure 8 
suggests that the lattice data are consistent with 
the ones from perturbation theory, although the 
simulation data have the tendency to be below 
the curves from perturbation theory. The same 

3 At Mh = 35 GeV only the value of Tc/Mh from lattice 
perturbation theory alone is given, assuming the same size 
of the error as at the other Higgs mass values 



is true for the situation in the 3-dimensional re- 
duced model, see table 1. 




Mh [GeV] 

Figure 5. Tc/Mh as a function of Mh ■ The solid 
lines are from continuum 2-loop computations [5]. 
The full symbols are from the 4-dimensional sim- 
ulation with square denoting Li = 3 and triangles 
Li = 2 results. The full circles are results from 
lattice perturbation theory extrapolated to the 
continuum and the error taken from the simula- 
tion data. 

In fig. 6 the situation for the important quan- 
tity vt/Tc is given. The scalar expectation value 
vt is in each case computed from the jump 
A of the gauge invariant quantity p"^, vt = 
^ A < > . Here lattice perturbation theory 
does not describe the data as well as for Tc/Mh 
(see fig. 4) and is therefore left out. In fig. 6 
results from simulations in the 3-dimensional re- 
duced model are included, denoted by the starred 
symbols. It can be seen that the different lat- 
tice simulation methods are compatible. In addi- 
tion, although the perturbatively obtained curves 
do not coincide with the Monte Carlo data, the 
agreement is surprisingly good, taken the general 
believe that perturbation theory might be inad- 
equate to describe the electroweak phase transi- 
tion. 
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20 40 60 80 

Mh [GeV] 

Figure 6. Same as fig. 5 for vt/Tc. In this fig- 
ure the results from the simulations in the 3- 
dimensional reduced model are included, denoted 
by starred symbols. 

6. Some open questions 

6.1. symmetric phase 

Although the results as summarized in the pre- 
vious section suggests that perturbation theory 
describes the electroweak phase transition reason- 
ably well, the situation for the physics of the high 
temperature symmetry restored phase is not as 
clear. In [55] an attempt was made to compute 
the spectrum in the symmetric phase within the 
3-dimensional SU(2)-IIiggs model. Starting with 
gap-equations, the phase transition itself and the 
behaviour of several observables close to the tran- 
sition could be described self-consistently. In par- 
ticular, it was found that in the symmetric phase, 
the effective model describing the system is again 
a SU(2)-IIiggs model with the form of eq.(3). 
Only the values of the couplings are different in 
the symmetric and the broken phases. This ap- 
pealing picture suggests that also the symmet- 
ric region can be treated perturbatively. If both 
regions are really analytically connected as indi- 
cated in fig. 1, the above scenario appears to be 
very plausible. 

The results from the work in [55] can be tested 



partly by Monte Carlo simulations. For exam- 
ple, the spectrum in the broken phase at finite 
temperature agrees very well between perturba- 
tion theory and numerical simulations. However, 
the situation in the symmetric phase is drasti- 
cally different. Here the quantities of interest are 
masses compared to the confinement scale g'^T . 
In [55] an upper bound for mw 1 9^T with mw 
the vector boson mass was given. The vector bo- 
son mass can also be determined in numerical 
simulations from suitable correlation functions. 
The results are shown in fig. 7. Here the values 
from simulations in 4 dimensions [40] and a re- 
cent simulation in the 3-dimensional Higgs model 
at Mh = 35GeV [56] are shown. The shaded area 
gives the allowed values from [55] with the upper 
bound for mw 1 9^T < 0.29. A clear discrepancy 
is encountered. Simulation results in the effective 
3-dimensional reduced model [27] are compatible 
with the numbers from the other two simulations 
as shown in the plot. 

At the moment it is not clear what the rea- 
son for this large difference is. The problem in 
the numerical simulations is that the correlation 
functions are very noisy in the symmetric phase 
and one might easily miss a low lying state as it 
disappears in the noise of the tail of the corre- 
lation functions. On the other hand, the picture 
as developed and described above in perturbation 
theory may not be valid. Clearly, more analytical 
and numerical work has to be done to clarify this 
important question. 

6.2. sphaleron rates 

For an estimate of the rate of baryon num- 
ber violating processes in the minimal standard 
model a knowledge of the sphaleron transition 
rate is very important. An earlier review of the 
work done on the sphaleron transition rate with 
emphasis on numerical simulations and with a 
much more comprehensive list of references can 
be found in [12]. In the low temperature phase 
the sphaleron rate describes how fast the baryon 
asymmetry is washed out after completion of 
the phase transition. Existing results for the 
sphaleron rate come from the semi-classical ap- 
proximation [10] and a non-perturbative test of 
these results is highly desirable. For this, how- 
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20 30 40 50 60 
Mh [GeV] 

Figure 7. The finite temperature vector boson 
mass divided by the confinement scale g^T. The 
shaded area indicates the allowed value as given 
in [55]. The simulation points are taken from [40] 
and from [56]. 

ever, one would need a real time simulation of 
the full quantum system. Unfortunately, such a 
method, up to now, does not exist. The best that 
can be done are real time simulations of the clas- 
sical system using classical Hamilton equations of 
motion. 

These classical real time simulations are per- 
formed by taking a Langevin type of algorithm 
to produce thermalized configurations [57]. One 
may alternatively use realistic heat bath meth- 
ods [59]. The thermalized configurations are then 
taken as input for the integration of classical 
Hamilton equations of motion. It is very im- 
portant that in the simulations the Gauss con- 
straints are respected. For the 2-dimensional 
abelian Higgs model, which has been investigated 
mostly so far, this is easy to implement. However, 
it was only recently that a method was found 
which respects the Gauss constraints exactly also 
for non-abelian gauge theories with or without 
matter fields [60]. 

The sphaleron rate F is related to the change 
A of the Chern-Simons number C during the real 
time simulation, A(t) = Tt where t denotes real 



time and 

A(t) = ^ I VpVqAC(t)e-"'^P-''^^'^ (23) 
with 

ACit) = [Cipit), qit)) - dpiO), qiO))f (24) 

is the canonical ensemble average at real time t 
with respect to the classical ensemble defined by 
the Hamiltonian H . 

For the low temperature regime of the 2- 
dimensional U(l)-Higgs model the semi-classical 
result has been worked out exactly [58]. The data 
from numerical simulations agree completely with 
the semi-classical ones. This is illustrated in fig. 8. 
The figure shows In F with F = F/m|_L and m$ 
the scalar mass in the U(l)-Higgs model. InF 
is plotted against the rescaled inverse tempera- 
ture l3' (see [57]). The figure suggests that at 
least for low temperatures the classical descrip- 
tion of the system is valid. This gives confidence 
in the semi-classical estimate for the transition 
rate shortly after the phase transition was com- 
pleted. It would be nice to see that the same is 
true for the 4-dimensional SU(2)-Higgs model and 
therefore to confirm eq.(l) which leads to relation 

At high temperatures where the system is in 
the symmetric phase, no exact analytical results 
are available. In fig. 8 the numerically obtained 
results are fitted to a behaviour. The data 
seem to respect this form, see however [59]. The 
two different curves correspond to two values of 
the lattice spacing indicating that the rate is lat- 
tice spacing dependent. The results as obtained 
in the 2-dimensional model are very promising. It 
seems to be possible in the near future that also 
classical simulations in the 4-dimensional SU(2)- 
Higgs model can be performed with the help of 
the method suggested in [60]. In 4 dimensions the 
sphaleron transition rate is expected to show a 
law, F = k(awT)^. Recent results indicate that 
for the pure gauge theory in four dimensions the 
behaviour of the rate is realized with a clas- 
sical value of the coefficient k = 1.09(4) [61]. A 
recent perturbative estimate in the 3-dimensional 
SU(2)-Higgs model reveals, however, k 0.01 
[62]. It remains to be seen, whether the inclu- 
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Figure 8. The sphaleron rate as a function of 
the inverse scaled temperature [3' [57]. For low 
temperatures the solid line is the exact result [58]. 
The curves at high temperature are fits, assuming 
a l//3'^ behaviour. 

sion of matter fields can change the coefficient k 
that drastically. 

6.3. crossover scenario 

It has been suggested by various groups doing 
numerical simulations [16], perturbation theory 
[55] and average action approach [34] that for a 
large enough Higgs mass Mh ~ 100 GeV the first 
order electroweak phase transition changes into 
an analytical crossover behaviour. It would be 
very important if one could clarify this picture. 
One way would be to repeat the old simulations 
which were done at -for todays standards- rather 
small lattices. An alternative would be to find an 
analyticity proof similar to the Osterwalder-Seiler 
proof [63] for the zero temperature system. 

6.4. linked cluster expansion 

As already mentioned in the text, the dimen- 
sional reduction program, performed on the level 
of the action, has been criticized. It would be im- 
portant to clarify the role of the non-local terms 
appearing when the reduction is done at the 2- 
loop level. Approaches to shed light on this 



problem or circumventing it are the concept of 
block spin transformation [32] or the matching 
of Green's functions [26]. Another interesting 
method is a recently performed high order linked 
cluster -hopping parameter- expansion [64]. The 
computation was performed up to the 18th order 
in N-component scalar field theories. As a first 
output, the method allows for a precise determi- 
nation of critical exponents. 

As a result, it could be shown that the (j)"^ model 
with # = 4 is governed by the 3-dimensional crit- 
ical exponents at high temperatures. This is in 
agreement with earlier lattice simulation results 
[17,18] and analytical investigations [19]. In ad- 
dition, it was found that the renormalized quartic 
coupling is weak close to the critical temperature 
and that (f)^' terms might be important. 

7. Conclusions 

In the introduction a moderate aim was formu- 
lated, namely, to confront results as obtained in 
perturbation theory with data from non-pertur- 
bative lattice simulations for Higgs masses Mh 
^ 70GeV. This goal has certainly been achieved 
by both, the 4-dimensional and the 3-dimensional 
reduced approaches. These two methods of in- 
vestigating the electroweak phase transition by 
means of numerical computations give compat- 
ible results providing an important cross check 
of the numerical data. The most surprising out- 
come from the numerical simulations is that up 
to Higgs masses of about 70GeV the phase transi- 
tion is well described by 2-loop perturbation the- 
ory. This is in contrast to earlier expectations 
that perturbation theory will not be able at all to 
say anything about the transition when physics of 
the symmetric phase has to be taken into account. 
It is a safe conclusion that we have a qualitative 
and quantitative good understanding of the finite 
temperature electroweak phase transition up to 
Higgs masses of about 70GeV. 

It is now well established that the electroweak 
phase transition is of first order if the Higgs mass 
is taken to be smaller than about 70GeV. The 
strength of the phase transition decreases rapidly 
with increasing Higgs mass. This is seen in ba- 
sically all physical observables, latent heat, sur- 
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face tension and jump of the scalar field expecta- 
tion value. The decrease is so pronounced that it 
becomes hard to imagine that for realistic Higgs 
mass values which are not ruled out by experi- 
ment the electroweak phase transition can pro- 
duce the right amount of baryon asymmetry. In 
order to complete the picture, simulations at a 
Higgs mass value of about the physical W-boson 
mass are mandatory. However, there the phase 
transition might be so weak that it will become 
very difficult to resolve its order [27]. 

Although, as mentioned above, the numerically 
obtained results are in reasonable agreement with 
the ones obtained in perturbation theory, there 
are slight discrepancies. In these differences the 
non-perturbative effects are presumably hidden. 
It seems that for an accurate quantification of 
the non-perturbative effects very precise numbers 
from numerical simulations are necessary. For 
this purpose computations in the 3-dimensional 
reduced model are probably the most appropri- 
ate. 

Open questions concern the high temperature 
regime. It is still not clear whether the phase 
transition will turn into the anticipated analyti- 
cal crossover behaviour at Higgs masses at about 
lOOGeV. The physics of the symmetric phase has 
to be clarified. This includes on the one hand the 
spectrum where a clear discrepancy between lat- 
tice simulations and perturbative calculations are 
encountered. On the other hand there is still a lot 
of work left to obtain a better understanding of 
the sphaleron rate. 
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